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Measurements of CKM angle 0i with charmless penguins at Belle 
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We present measurements of time-dependent CP violation parameters in B° — > K°Tr°, B° —* K^n^ir" and B° —* 
KgK + K~ decays. The latter two are extracted using time-dependent Dalitz plot analyses. These results are obtained 
from a large data sample that contains 657 X 10 6 BB pairs collected at the T(4S) resonance with the Belle detector 
at the KEKB asymmetric-energy e + e _ collider. 
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Q 1. INTRODUCTION 

O ■ 
(N ' 

CP violation in the Standard Model (SM) arises from an irreducible complex phase in the Cabibbo-Kobayashi- 
O ■ Maskawa (CKM) quark-mixing matrix [l|, |2j. Of recent interest is CP violation in b — > qqs transitions which 
proceeds by loop diagrams that may be affected by new particles in various extensions of the SM. Furthermore, the 
' CP asymmetries in b — > qqs transitions are predicted in the SM to be slightly higher than those observed in b — > ccs 
transitions. However, current experimental measurements |3j tend to be lower than those for b — > ccs transitions 
r— i motivating more precise experimental determinations. 

The decay of the T(4S) produces a BB pair of which one (Br cc ) may be fully reconstructed while the other (-Bxag) 
may reveal its flavour. The proper time interval between B^ cc and B^ ag is defined as At = tjx cc — t^ ag and from 
coherent BB production in the T(45) decay, the time-dependent decay rate for a quasi- two-body mode when £?xag 
■ possesses flavour q (B°: q = +1, B°: q = — 1), is given by 
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' where Tb° is the B° lifetime and Arrid is the B°B° mass difference. This assumes no CP violation in mixing, \q/p\ = 1, 
and that the B a B° lifetime difference is negligible. The parameter, Acp, denotes the direct CP violating component 
and Sep represents mixing-induced CP violation. For time-dependent Dalitz plot analyses, the time-dependent 

, decay rate is written as 
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The Dalitz-dependent amplitudes, A, can be written in the isobar approximation as a superposition of intermediate 
decay channels, i, 

A(s+,s_) = ^a l F l (s + ,s_), A(s-,s+) = ^a 4 ^(s_,s + ), (3) 

i i 

where = aj(l + Ci)e l ^ bi+di ^ for A and di = — Ci)e^ bi ^ di ^ for A are complex coefficients describing the rela- 
tive magnitudes and phases between the decay channels. The form factors, i^(s+,s_), depend on the Dalitz plot 
coordinates, s±, and describe the invariant mass and angular distribution probabilities. For a CP eigenstate, i, the 
time-dependent CP violation parameters can be calculated as 



i , I — 2c» , c ... _ — 23(a 4 a*) l — c i , ll 

A CP(V = !2 I HI 13 = and " W5CP(*) = i | 2 + |- |2 = YT^2 sm2 ^i ( J ) ( 4 ) 



and </>i ff (i) = arg(oiO*)/2 = di is directly accessible as a fit parameter. Consequently, both Acp{i) and Scp{i) are 
restricted to reside in the physical region. The relative fraction of each component can be calculated with, 

= {\a i \ 2 + \a i \ 2 )^F i { 3+ , S -)F:' S+ , S -)ds + d S - 
h j\\A\ 2 + \A\ 2 )ds+ds- ' 1 ' 
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2. DATASET, DETECTOR AND BASIC ANALYSIS TECHNIQUE 

These measurements of CP violating parameters are based on 657 x 10 6 BB pairs collected with the Belle detec- 
tor Q at the KEKB asymmetric-energy e + e _ (3.5 on 8 GeV) collider Operating with a peak luminosity that 
exceeds 1.7 x 10 34 cm _2 s -1 , the collider produces the T(4S) resonance (y/s = 10.58 GeV) with a Lorentz boost of 
/?7 = 0.425, opposite to the positron beam direction, z, which usually decays into a BB pair. 

Reconstructed B candidates are described with two kinematic variables: the beam-constrained mass, M\, c = 
V^bcam) 2 " bs MS ) 2 and the energy difference, AE = £§ MS - ££MS where e cms ig the beam energy and e cms 
(Ps ) is the energy (momentum) of the B meson all evaluated in the centre-of-mass system (CMS). The dominant 
background in the reconstruction of B^ cc is from continuum (e + e~ — > qq) events. Since their topology tends to be 
jet-like in contrast to the spherical BB decay, continuum can be suppressed with a Fisher discriminant based on 
modified Fox- Wolfram moments This discriminant is combined with the polar angle of the B candidate in the 
CMS to form a likelihood ratio, TZ-s/Bi which separates continuum from BB events. 

Since the B° and B° mesons are approximately at rest in the T(45) Enter-of-Mass System (CMS), the difference 
in decay time between the BB pair, At, can be determined from the displacement in z between the final state decay 
vertices, At ~ (zr cc — ^Tag) / 'ftjc = Az/fijc. To obtain the At distribution, we reconstruct the tag-side vertex from 
the tracks not used to reconstruct B^ ec Q and employ the flavour tagging routine described in Ref. Q . The tagging 
information is represented by two parameters, the -B^-ag flavour, q and r. The parameter, r, is an event-by-event, 
MC determined flavour-tagging dilution factor that ranges from r — for no flavour discrimination to r — 1 for 
unambiguous flavour assignment. 



3. TIME-DEPENDENT CP VIOLATION MEASUREMENT IN B° -> K°n° 

Direct CP violation has been observed in B° — * K + tt~ [h]] and is found to be significantly different from that 
in B^ — > K^ir® . This unexpected result may indicate the presence of new physics (NP) or poor understanding of 
strong interaction effects in B decays. A model-independent test for NP is possible via an isospin sum rule with high 
statistics [ll| which gives a relation among the direct CP violation asymmetries measured in all possible B — > Kn 
modes, i.e. B° -> K°ir°, K+ir- , B+ -> K+ir° and K°ir+, 

B(K°n+)r B o 2B(K+n°)T B a 2B(K°n°) 

A K+n - +AK°« +B{K+n _ )TB+ - a ^ b{k+7T - )tb+ + A ^" B{K+7r -y 

Here, B represents the branching ratio of each decay mode and t b + is the lifetime of the charged B meson. The sum 
rule's theoretical precision is determined by SU(2) flavour symmetry, i.e. a few %, therefore the sum rule provides 
a clean test for new physics. As a fundamental input, the CP violation measurement in B° — > K°ir° is currently 
the least known, experimentally. Since the branching fractions and CP asymmetries of other B — » Kir decay modes 
have been measured to good precision [3], Ak ^ is predicted with a small error. In addition to the B° — > Kgir 
mode, we measure CP asymmetry in B° — > K^tt decay for the first time, in order to maximize sensitivity to the 
direct CP violation parameter, Ak tt°- 

The signal yield of B° — > Kgir events is found from a three-dimensional extended unbinned maximum likelihood 
fit to Mbc, AE and 72-g/B to be 657 ± 37 events. For B° — > K l tt°, M^ c can be calculated from the direction of the 
Kl cluster while AE cannot be calculated. The signal yield of 285 ± 52 events is extracted from a two-dimensional 
fit to Mbc and TZs/b an d has a significance of 3.7<7 including systematic uncertainties. 

As the vertex position of B° — > K^ir cannot be determined and the vertex reconstruction efficiency of B° — > K < s 7r ° 
is ~ 33% due to the long lifetime of the Kg, these events can still be used to calculate Acp by integrating Eq. Q]over 
At. We extract the CP parameters, 

A C p = +0.14 ± 0.13 (stat) ±0.06 (syst), 

S C p = +0.67 ±0.31 (stat) ±0.08 (syst) (7) 
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Figure 1: The left plot shows the background subtracted fit to At for the good tags region, 0.5 < r < 1.0, where the solid 
(dashed) curve represents B° (B°) tags. The right plot shows the B°B° raw asymmetry, (N B o — N B o)/(N B o + N B o), where 
N B o (N B o) is the number of signal B° (B°) tags in At. 



and the fit results for the At component is shown in Fig.[TJ We find that the mixing induced component is consistent 
with charmonium, Scp(b — > ccs) — 0.681 ± 0.025 and that there is a 1.9a deviation between our measurement of 
Acp and the expectation from the isospin sum rule, Acp — — 0.17 ± 0.06 [11|. 



4. TIME-DEPENDENT DALITZ PLOT CP VIOLATION MEASUREMENT IN B c 



The signal yield of B° — ► Kgir + ir events is found from a one-dimensional fit to AE to be 1944 ± 98 events and 
the resonances considered in the signal model (Eq. EJ) are the ^±(892), ^±(1430), p°(770), /o(980), / 2 (1270), 
/x(1300) and a non-resonant component. We obtain the CP parameters, 



Solution 1: -21og£ = 18472.5 

f(K* + {U30)ir-) = 61.7 ± 10.4% 

A c p(p° (770) K° s ) = +0.03tg;g± 0.11 ±0.10, 

4f {p°(770)K° s ) = (+20.0l|j±3.2±3.5)°, 

A C p(fo(980)K^) = -0.06 ±0.17 ±0.07 ±0.09, 

<j>f(fo(980)K° s ) = C+12.7±l;| ±2.8 ±3.3)°, 



Solution 2: -21og£ = 18465.0 

/(A'* + (1430)7r-) = 17.4 ± 5.0% 

A C p{p°{77Q)K%) = -0.16 ±0.24 ±0.12 ±0.10, 

<j>f{p°{770)K G s ) = (+22.8 ±7.5 ±3.3 ±3.5)°, 

A C p(fo(980)K%) = +0.00 ±0.17 ±0.06 ±0.09, 

cj ) f(fo(980)K° s ) = (+14.8^; 3 7 ±2.7±3.3)°, (8) 



where the first error is statistical, the second is systematic and the third is the Dalitz plot signal model uncertainty. 
Figure [2] shows the fit results. The high -K~q + (1430)7t - fraction of Solution 1 is in agreement with some phenomeno- 
and may also be favoured by the total K — n S'-wave phase shift when compared with that 



logical estimates [1 
measured by LASS 

solutions where Acp is consistent with null asymmetry and cf>1 s agrees with charmonium. 



14j . As the likelihood difference is not found to be significant, we do not rule between these 



5. TIME-DEPENDENT DALITZ PLOT CP VIOLATION MEASUREMENT IN B° -> K S K + R- 

The signal yield of B° — > KgK + K~ is extracted to be 1269 ± 51 events from a two-dimensional fit to M^ c and AE 
in each r-bin. We consider the /o(980), 0(1020), /x(1500), Xco, and a non-resonant component in the signal model. 
Four solutions with similar likelihood are found where two solutions arise from the interference between /o(980) and 
the non-resonant component and the other two from /x(1500) and the non-resonant component. Using external 
information from B° -> i^Tr+Tr - , if the /x(1500) is the / (1500) for both B° -> K^n+n- and B° -> K%K + R-, 
the ratio of branching fractions, S(/o(1500) — > 7r + 7r _ )/S(/o(1500) — > K + K~), prefers the solution with the low 
/x(1500)^ fraction. Similarly the ratio, B(/ o (980) 7r+7r-)/(S(/ (980) -> tt+tt") + B(/ (980) K + K~)) 
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Figure 2: Time-dependent Dalitz plot fit results for B° — > K%-k + -k~ in (a), the p°(770) region and (b), the /o(980) region. 
The top plots show the At distribution for B° (solid) and B° (dashed) tags. These plots contain only good tags, 0.5 < 
r < 1.0 and the dotted curve represents the background contribution. The bottom plots show the B°B° raw asymmetry, 
(iV s o - N S o)/(N B o + Ng ), where JV B o (N S o) is the number of B° (B°) tags in At. 



prefers the solution with the low fo(980)Kg fraction. The preferred solution is, 

A C p(fo{9S0)K° s ) = -0.02 ± 0.34 ± 0.08 ± 0.09 

<Pf(fo(?80)K%) = (28.2±|J ± 2.0 ± 2.0)° 
A C p (</>(1020)if§) = +0Mto 2 2 l ± 0.04 ± 0.09 

<t>f {(j)(1020)K o s ) = (21.2^ 8 4 ± 2.0 ± 2.0)° (9) 

where the first error is statistical, the second is systematic and the third is the Dalitz plot signal model uncertainty. 
Figure [3] shows the fit results where Acp is consistent with null asymmetry and c^ ff agrees with charmonium. 
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Figure 3: Time-dependent Dalitz plot fit results for B° -> K%K + K~ in the 0(1020) region. The left plot show the At 
distribution for B° (blue) and B° (red) tags for good tags, 0.5 < r < 1.0. The right plot shows the B°B° raw asymmetry 
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